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Abstract. We use J, H and Ks 2MASS photometry to study colour-magnitude (CMDs) and colour-colour diagrams, 
£ — 1 structure and mass distribution in the ionizing open cluster NGC 6611. Reddening variation throughout the cluster 

region is taken into account followed by field-star decontamination of the CMDs. Decontamination is also applied 
to derive the density profile and luminosity functions in the core, halo and overall (whole cluster) regions. The 
field-star decontamination showed that the lower limit of the main sequence (MS) occurs at ~ 5 Mq . Based 
on the fraction of Ks excess stars in the colour-colour diagram we estimate an age of 1.3 ± 0.3 Myr which is 
consistent with the presence of a large number of pre-main sequence (PMS) stars. The distance from the Sun was 
estimated from known O V stars in the cluster area and the turn-on stars connecting the PMS and MS, resulting 
in d0 = 1.8±0.5kpc. The radial density distribution including MS and PMS stars is fitted by a King profile with 
a core radius R CO ro = 0.70 ± 0.08 pc. The cluster density profile merges into the background at a limiting radius 
ly-^ ' Riim = 6.5 ± 0.5 pc. From the field-star subtracted luminosity functions we derive the mass functions (MFs) in 

the form <j>(m) oc m ( -~ 1+x \ In the halo and through the whole cluster the MFs have slopes x = 1-52 ± 0.13 and 
X = 1.45±0.12, respectively, thus slightly steeper than Salpeter's IMF. In the core the MF is flat, \ — 0.62 ±0.16, 
indicating some degree of mass segregation since the cluster age is a factor ~ 2 larger than the relaxation time. 
Because of the very young age of NGC 6611, part of this effect appears to be related to the molecular cloud- 
fragmentation process itself. We detect 362 ± 120 PMS stars. The total observed mass including detected MS (in 
the range 5 — 85 Mq) and PMS stars amounts to ~ I6OOM0, thus more massive than the Trapezium cluster. 
Compared to older open clusters of different masses, the overall NGC 6611 fits in the relations involving structural 
and dynamical parameters. However, the core is atypical in the sense that it looks like an old/dynamically evolved 
core. Again, part of this effect must be linked to formation processes. 
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1. Introduction stellar-mass distribution in young star clusters may shed 

light on this issue. A suitable target for this kind of study 

The distribution of stellar masses at star-cluster birth is is the young open cluster NGC 6611, which is responsi- 

a complex process not yet fully understood. Recent ob- ble fo r the ionization of th e Hll region Sh2-49 (Sharpless 

servations of the stellar content of star-forming regions flT^ . aka Gum 83 (Gum 115551 or RCW165 (Rodgers, 

in molecular clouds, rich star clusters and the Galactic Campbell, & Whi teoak ITrMTI) as seen in the optical, or 

field suggest that the initial mass function (IMF) has W37 (Westerhout ITrj53|) in the radio, 
similar properties in these very different environments 

(Kroupa 2002). This seems to indicate that the initial The open c i uste r NGC 6611 (M 16, Mel-198, Cr375, 

distribution of stellar masses should depend only on the OC154 - Alter, Ruprecht, & Vanisek 1117701 consists of a 

process of molecular cloud fragmentation, and that the syste m of early-type stars embedded in emission nebulos- 

fragmentation would have to produce similar IMFs de- ity ( The Eagle Nebula) with spatially variable dust atten- 

spite very different initial conditions, a physic al pro cess uation in the op tical range. Structurally the region shows 

which still lacks a better description (Kroupa |H2). In some ne bulosity features, the so-called elephant trunks 

this sense, detailed analysis of the spatial structure and ( the "Pillars of Creation", Hester, Scowen, Sankrit, et al. 

at the tips of which new-born stars start to be vis- 
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The cluster is located at I rj 17° and b » 0.8°, in 
the Sagittarius spiral arm (Georgelin, & Georgelin^)70), 
a region presenting conspicuous dense dust clouds. The 
whole region is an active site where star-formation pro- 
cesses seem to occur giving birth to low to high-mass stars 
(Kroupa 2004). Evidence has been presented in recent 
years that NGC6611 shows an extremely flat MF (e.g. 
Sagar & Joshi lT^I Sagar et al. 115551 Hillenbrand et al. 
H333I Massey et al.UMB- The distance of NGC 6611 is de- 
pendent on the extinction law but several works constrain 
it to d = 2.3 ± 0.3 kpc (The et al. IT33U1 Hillenbrand 
et al. ITMSl Belikov et al lTMjjl . The cluster stars have 
been observed in evolutionary stages which are not recon- 
cilable with a single age system (de Winter et al. 1997). 
Also, isochrone age determinations are ambiguous for such 
young clusters because of the lack of constraints in the 
CMD morphology. However, several studies show that the 
cluster is certainly younger than lOMyr with the youngest 
members still being formed (e.g. Hillenbrand et al. [1993 
de Winter HM7I Belikov et alEOOO). NGC 6611 presents a 
significant number of PMS stars, de Winter et al. (1997) 
revised the subject of PMS membership and related issues 
by studying PMS candidates individually via photometry 
and spectroscopy. They find that the sample objects are 
not coeval, with an age spread of w 6Myr being necessary 
to explain the cluster stellar content. 

Our main goals with the present work are (i) carry out 
a detailed analysis of the spatial structure of NGC 6611, 
(ii) analyze the spatial variation of the cluster MFs, (in) 
derive properties of the MS and PMS stars, and (iv) infer 
on the dynamical state of this very young open cluster. 
For the sake of spatial and photometric uniformity, we 
employ J, H and K s 2MASS 1 photometry. The 2MASS 
Point Source Catalogue (PSC) is uniform reaching rela- 
tively faint magnitudes covering nearly all the sky, allow- 
ing a proper background definition even for clusters with 
large angular sizes (e.g. Bonatto, Bica, & Santos Jr. 2004 
Bonatto, Bica, & Pavani l2HU3jl . 

This paper is organized as follows. In Sect. |3 we con- 
dense previous results on NGC 6611. In Sect. we present 
the 2MASS data, correct the photometry for differential 
reddening, subtract the field-star contamination and ana- 
lyze the CMDs and colour-colour diagrams. In Sect.0]we 
discuss the cluster structure. In Sect. |5] we derive LFs and 
MFs and discuss stellar content properties. In Sect. El we 
compare NGC 6611 with nearby, older open clusters and 
discuss dynamical states. Finally, concluding remarks are 
given in Sect. 

2. The young open cluster NGC 6611 

The equatorial young open cluster NGC 6611 is included 
in the Ser OB 1 association and the W37 molecular cloud. 
It is part of a larger sequence of star-forming events which 

1 The Two Micron All Sky Survey, All Sky 
data release (Skrutskie et al. I1997II . available at 
://www. ipac.caltech.edu/2mass/releases/allsky/ 



have just begun in the region and are still embedded in 
the molecular cloud (Hillenbrand et al. 11993(1 . The cluster 
contains 8 OV-stars (Bosch, Morrell, & Niemela, 1999). 

The first comprehensive photometric study on this ob- 
ject was carried out by Walker tl96T|l who obtained UBV 
photoelectric and photographic observations of NGC 6611 
limited at V=16.7. He found that in the southern part of 
the region the colour excess is relatively small, E(B — V) = 
0.6, while to the north the reddening increases irregularly 
reaching E(B — V) = 1.73. Stars above the MS were con- 
sidered possible members in the gravitational contraction 
phase. He derived an age of 1.8 Myr based on the analysis 
of 532 probable member stars. 

Sagar & Joshi (1979) derived a distance to the Sun 
d Q = 3.2 ± 0.3 kpc, an age in the range 3.0 - 5.5 Myr, 
and observed an extremely flat luminosity function in this 
open cluster. 

Based on photoelectric UBV photometry, Sagar et al. 
( 198<3) derive a flat MF with a slope \ ~ —0.2 for stars in 
the mass range 9 < m(M ) < 70. 

By means of photometric observations over the range 
from 0.3 to 3.7 ^m Chini & Wargau (J 990) provide evi- 
dence of an abnormal extinction law only at wavelengths 
shorter than 0.55 /im. They conclude that all PMS stars 
previously identified by Walker 1(1961(1 present reddening 
due to normal dust, therefore suggesting that they are 
foreground objects. 

Hillenbrand et al. l(TM3|l combined CCD UBV obser- 
vations with near-infrared JHK ones to establish a theo- 
retical HR diagram for the cluster, from which they con- 
clude that NGC 6611 is actively forming 3 - 8M© stars. 
Using the optical data they infer from a spatial coverage 
of w 1 600' 2 an overall MS MF slope of x = 0.1 - 0.3 in 
the mass range 5 < m(M©) < 85. It is worth noting that 
only 4 stars have masses between 40 and 85 M© , for which 
they derive an age of 2±1 Myr. They also present evidence 
that one 30 M© star has been formed ~ 6 Myr ago. They 
derive d Q = 2.0 ±0.1 kpc. 

By constructing new HR diagrams from the data pub- 
lished by Hillenbrand et al. 1)1993(1 and applying a different 
mass calibration, Massey et al. 1(19951 derived an age in 
the range 1-5 Myr, \ — —0.3 ± 0.2, and a stellar density 
of p = 5.4 x 10 4 kpc~ 2 for stars above 10 M©. 

Belikov et al. 1(1999(1 analyzed the spatial and proper 
motion distribution of cluster stars, assuming a normal 
distribution for the spatial velocities. They derived a dis- 
tance to the Sun d© = 2.14±0.10 kpc and obtained spatial 
dispersions of off 5 = 1.20' and <j™™ na = 4.77', respec- 
tively for the core and corona subsystems in NGC 6611. 
They used the same data to measure the differential red- 
dening throughout the cluster area. 

Belikov et al. l(2DUU)> studied NGC 6611 based on a 
compilation of photographic, photoelectric and CCD op- 
tical (UBV) observations covering a circular area of 45.6' 
in diameter centered on NGC 6611. They present radial 
profiles of the star density distribution in different az- 
imuthal sectors. Small differences found between the den- 
sity of the north-west and south-east sectors were at- 
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Fig. 1. Left panel: XDSS R image. Right panel: 2MASS J image. The area covered in both images is 1° x 1°. The 
small circle at the center of each image indicates the core, and the large one indicates the size of NGC 6611 (limiting 
radius). 



tributed to irregularities of the absorption due to inter- 
stellar dust rather than to stellar density fluctuations. The 
south-east sector of the cluster is less affected by dust 
than the north-west one. The total radial density distri- 
bution allowed them to conclude that virtually no clus- 
ter members are present beyond the edge of the cluster 
corona, at 3a^Y° na = 14.3'. They derive an age of 6Myr 
for NGC 6611. Based on the number frequency of member 
stars in the range 2.1 < m(M Q ) < 85 they infer a flat MF 
slope, x — 0-2. 

Tadross et al. (2002) presented an analysis of UBV 
CCD observations intending to derive morphological pa- 
rameters for a large sample of open clusters. For NGC 6611 
they obtained a colour excess E(B — V) = 0.68, d© « 
1.6 kpc, age w 3.2 Myr, linear diameter D = 3.6 pc, 
Galactocentric distance dec = 6.94 kpc, number of mem- 
ber stars N* = 282, and a total mass of hims+pms = 
435 M Q . 

The presence of significant numbers of PMS stars in 
NGC 6611 has been confirmed by means of the observa- 
tion of a large number of emission-line stars and IR disk 
signatures (Hillenbrand et al. 119931 Massey et al. Q.995), 
which is an expected population in such a young cluster. 
Recently Oliveira et al. lf2T705 ) obtained IZJHKL' photom- 
etry aiming at identifying PMS stars in NGC 6611. Their 
Fig. 2 shows the CMDs I x (I — Z) for the central area 
(7' x 6') and for a control field 16' away from the cluster 
center. The PMS population, defined within 12 < I < 19, 
corresponding to the mass range 0.4 < m(M Q ) < 7, stands 
out in the central area and decreases significantly towards 
the control field. The MS population does not show such 



trend, its spatial distribution remaining apparently con- 
stant throughout the fields. 

In the left panel of Fig.""f]we provide an R XDSS image 
encompassing an area of 1° x 1° centered on NGC 6611, 
extracted from the Canadian Astronomy Data Centre 
(CADC 2 ). In the right panel we show the corresponding 
2MASS J image. The stark contrast between the R and J 
images reflects primarily the effect of the Ha emission in 
the R band and how dust is much more penetrable in the 
near-infrared bands. Thus, the analysis of NGC 6611 and 
its outskirts with near-infrared data is naturally justified 
because it allows access to low-mass stars with moderate 
aperture telescopes. In this sense, Hillenbrand et al |1993 ) 
investigated NGC 6611 by means of deep near-infrared ob- 
servations covering a 15' x 15' area with useful limits of 
J=16, H=14.5 and K=13.5, similar to those of 2MASS. 

According to the WEBDA 3 open cluster database 
fMermilliod 11996*)! . the central coordinates of NGC 6611 
are (J2000) a = 18 ,l 18 m 48 s , and 5 = -13°48'24". 
However, the corresponding radial density profile (Sect.0} 
presented a dip for R = 0'. Accordingly, we searched 
for a new center by examining histograms for the num- 
ber of stars in 1' bins of right ascension and declina- 
tion. The resulting coordinates which maximize the den- 
sity of stars at the center are (J2000) a = 18M8 m 40.8 s , 
and 5 = -13°47'24.0", corresponding to I = 16.95° and 
b = 0.83°. In what follows we refer to these optimized 
coordinates as the center of NGC 6611. WEBDA gives a 



2 http://cadcwww.dao.nrc.ca/ 

3 http:/ / obswww. unige. ch/webc 
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Fig. 2. Spatial distribution of the number of stars 
(brighter than J = 15.8) around the central position of 
NGC 6611. The density of stars increases to the East and 
South, lcr Poisson errors are shown. 



colour excess of E(B — V) = 0.768, distance to the Sun 
d Q = 1.75 kpc, and age ss 7.6 Myr. 

3. Photometric parameters from 2MASS data 

The VizieR 4 tool was used to extract J, H and K s 2MASS 
photometry of the stars present in a circular area with 
radius R = 60' centered on the optimized coordinates de- 
scribed in the previous section. The faint-magnitude limit 
of the extracted stars is brigther than that correspond- 
ing to the 99.9% Point Source Catalogue Completeness 
Limit 5 , J = 15.8, H = 15.1 and Ks = 14.3, respectively. 

As a first step in the investigation of NGC 6611 we use 
the 2MASS data to examine in Fig. [5] the distribution of 
stars in the area surrounding the cluster. Each panel de- 
picts the number of stars brighter than J = 15.8 in strips 
of 5' width in declination (left panels) and right ascension 
(right panels) with respect to the optimized cluster center. 
The rather concentrated structure of NGC 6611 is appar- 
ent in panels (e) and (f). Within uncertainties the cluster 
appears to be spatially asymmetric, being more extended 



4 http://vizier.u-strasbg.fr/viz-bin/VizieR?-source= 11/246 

5 According to the Level 1 Requirement, according to 

http://www.ipac. caltech. edu/ '2mass /releases/ allsky /doc/ sec6=5al .html 
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Fig. 3. J x (J-H) CMDs of the central 9' of 
NGC 6611. Left panels: observed photometry. Middle pan- 
els: reddening-corrected CMDs according to Belikov et al. 
(1999). Right panels: reddening-corrected and field-star 
decontaminated CMDs. Note the extended MS and con- 
spicuous PMS population, which is restricted essentially 
to the central 6'. Field stars are shown as x in the middle 
panels. Filled circles in the bottom panels: O V stars from 
Bosch, Morrell & Nicmclii 



in right ascension (apparent diameter D a ~ 21', (panel 
(e)) than in declination (D,5 ~ 12', (panel (f)). The spatial 
distribution of stars around the cluster is non-uniform, in- 
creasing significantly above the lcr Poisson errors towards 
the East (panels (a), (c), (e) and (g)) and less markedly to 
the South (panels (b) and (d)). Part of this non- uniformity 
can be accounted for by differential reddening, which is im- 
portant throughout the cluster area (Belikov et al.[l999). 

3.1. Spatial variation of reddening 

Because of the low Galactic latitude (b — 0.83°), the field 
of NGC 6611 presents significant field-star contamination. 
This can be seen in the left panels of Fig.[3]where we show 
the observed J x (J — H) CMDs of three consecutive cir- 
cular regions (AR = 3') centered on NGC 6611. Besides 
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conspicuous MS stars, the integrated CMD up to R = 9' 
(bottom panel) presents important field-star contamina- 
tion. 

As discussed in Belikov et al. (1999}, the spatial field 
of NGC 6611 is critically affected by differential reddening 
with important variations both in E(B — V) and Ry in ar- 
eas separated by distances as small as ~ 2'. They present 
E(B — V) and Ry values derived from spectroscopy and 
proper motion of cluster stars in 2' x 2' cells covering an 
area of about 40' x 36' centered on NGC 6611. Note that 
E(B — V) in each cell corresponds to the absolute redden- 
ing value. 

We applied the above procedure to produce reddening- 
corrected 2MASS photometry for all stars in each cell to 
yield a total region reaching R w 23' from the center of 
NGC 6611. In the vacant cells we followed their suggestion 
and adopted the average values of E(B — V) and Ry- For 
reddening and absorption transformations in the 2MASS 
bands we use the relations Aj = 0.276 x R V E(B - V), 
A H = Aj-0.33E(B- V) and A Ks = Aj -0.49E(B - V), 
according to Dutra, Santiago & Bica ( 2002 , and references 
therein). 

The resulting reddening-corrected CMDs are shown in 
the middle panels of Fig. [3] Comparison of the reddening- 
corrected CMDs with the observed ones shows that 
the average colour-excess throughout the cluster area is 
(E(J - H)) = 0.25, corresponding to (E(B - V)} = 0.8. 

3.2. Field-star decontamination 

Finally, to uncover the intrinsic CMD morphology we ap- 
ply a field-star decontamination procedure. Because of 
the spatial limitation on the useful area caused by the 
differential-reddening correction (Sect. 13. If) , we consider 
as offset field the ring located at 17' < R < 22'. This area 
is large enough to produce statistical representativity of 
the field stars, both in magnitude and colours. 

To illustrate the amount and colour-magnitude distri- 
bution of the contamination in the area of NGC 6611 we 
plot in the middle panels of Fig.[3]the corresponding (same 
area) field-star contributions taken from rings with exter- 
nal radius at 22'. As a consequence we can see that most 
of the faint stars - in particular those mimicking a low-MS 
- are in fact field-stars. 

Based on the spatial number-density of stars in the 
offset field, the decontamination procedure estimates the 
number of field stars which within the ler Poisson fluc- 
tuation should be present in the cluster field. The ob- 
served CMD is then divided in colour/magnitude cells 
from which stars are randomly subtracted in a number 
consistent with the expected number of field stars in the 
same cell. The dimensions of the colour/magnitude cells 
can be subsequently changed so that the total number of 
stars subtracted throughout the wole cluster area matches 
the expected one, within the lcr Poisson fluctuation. Since 
the field stars are taken from an external annulus of fixed 
dimensions, corrections are made for the different solid 
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Fig. 4. (H - Ks)x(J - H) diagram of the central 3' (Panel 
(a)) and 9' (Panel (b)) of NGC 6611. Dashed line: 4Myr 
Padova isochrone. Dotted line: 1 Myr sequence of PMS 
stars (Siess, Dufour, & Forestini l2UUUjl . Solid line: MS 
stars from Schmidt-Kaler l|1982f) . Filled triangles: OV 
stars (Bosch, Morrell, & Niemela 1999). Filled circles: red- 
dening vectors; AE(B — V) = 0.25 between 2 consecutive 
circles. The corresponding field stars are shown in Panel 
(a) as 'x'. Both diagrams were produced after field-star 
decontamination. 

angles of cluster and offset field. This procedure can be 
applied to the cluster region as a whole or internal regions 
as well. Note that because this procedure actually excludes 
stars from the original files - thus artificially changing both 
the radial distribution of stars and the LF - we use it only 
for the sake of uncovering the intrinsic CMD morphology. 

The resulting field-star decontaminated (differential- 
reddening corrected) CMDs are shown in the right panels 
of Fig. |21 As expected, most of the faint-star distribution 
was subtracted by the decontamination procedure. The 
central region (0' < R < 3') presents an extended, vertical 
MS reaching from J as 11.5 to J w 5.5 which reflects the 
young age of NGC 6611. This region also presents a con- 
spicuous population of PMS stars (0.54 < (J — H) < 1.1) 
which extends essentially up to R ss 6'. The entry point 
(turn-on stars) of the PMS into the MS occurs at J « 11.1. 
The region 6' < R < 9' is sparsely populated by cluster 
stars, with a small fraction of MS and PMS stars still dis- 
cernible. 
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3.3. Age 

Deriving a precise age for clusters younger than ~ 10 Myr 
by means of isochrone fit is a difficult task because of 
the nearly-vertical MS and lack of post-MS stars. As a 
consequence of the scarcity of observational constraints in 
the CMD, several age-solutions are possible. In the case 
of NGC6611, the CMD (Fig.0 can be acceptably fitted 
with isochrones with ages up to 10 Myr. In the present case 
we use the set of solar-metallicity isochrones of Padova 
(Girardi et al. IMI^ computed with the 2MASS J, H 
and Ks filters 6 . The 2MASS transmission filters produced 
isochrones very similar to the Johnson ones, with differ- 
ences of at most 0.01 in (J — H) (Bonatto, Bica & Girardi 

mm. 

On the other hand, properties of young/embedded ob- 
jects, in particular the age, can be alternatively inferred 
by means of infrared colour-colour diagrams (e.g. Lada, 
Alves & Lada l99.6|. The fraction of stars with Kg excess 
in embedded clusters is an age indicator (e.g. Lada, Alves 
& Lada ll99rJl Soares & Bica l2002|) in the sense that larger 
fractions correspond to younger evolutionary stages. In 
Fig. 0] we show the (H — Ks) x ( J — H) diagram of the 
central 3' (Panel (a)) and 9' (Panel (b)) of NGC6611. 
Both diagrams were produced after differential-reddening 
correction and field-star decontamination. To compare 
stars in different evolutionary stages we show in Panel (a) 
the loci occupied by MS stars according to the youngest 
Padova isochrone (4 Myr) and to those of Schmidt-Kalcr 
((1112). The 1 Myr sequence of PMS stars (Siess, Dufour, & 
Forestini 2000 ) , and the loci of the O V stars with available 
spectroscopic data of Bosch, Morrell, & Niemela (|l999) 
are shown as well. Extinction effects can be estimated by 
the reddening vectors. 

For NGC6611 we measured within R < 3' a fraction 
of f cxc = 0.54±0.05 of Ks excess stars. Considering uncer- 
tainties this value is intermediate between those of the p 
Ophiuchi cluster (f exc = 0.5 - 0.7 - Greene & Young lTM21 
Strom, Kepner, & Strom I1995fl with an age of 1 Myr or 
less (Greene & Mever IT995|I and that of the Taurus Dark 
Nebula group (f cxc = 0.5 - Kenyon & Hartmann [1995) 
with an age of 1.5 Myr. Based on this we estimate an age 
of 1.3 ± 0.3Myr for NGC 6611. These values are similar to 
those of Trapezium, NGC 2327 cluster and BRC 27 cluster 
(Soares & Bica lMi^jl . 

This estimate is consistent with (i) the presence 
of a significant population of PMS (e.g. Walker 119611 
Hillenbrand et al. HMfl Massey et al. ITM>1 Oliveira et 
al. 11003), (H) the fraction of - 58% of PMS having cir- 
cumstellar discs (e.g. Hillenbrand et al. 119931 Massey et 
al. 119951 Oliveira et al. I2005J) . and (Hi) the fact that 
NGC 6611 is ionizing the Hil region Sh2-49 (e.g. Bosch, 
Morrell, & Niemela ITM?1 Hillenbrand et al. H^5|l . 



http:/ /pleiadi.pd. astro .it/isoc_photsy s . 1 /isoc_photsy s . 1 . html 



3.4. Distance from the Sun 

Similarly to the age, the lack of evolved features in the 
CMD of NGC 6611 precludes a precise distance determi- 
nation based directly on isochrone fit. In this sense we 
will derive the distance of NGC 6611 from the Sun us- 
ing two independent methods. Bosch, Morrell & Niemela 
( 1999J derived spectroscopic data of 8 MS O member stars 
(which are identified in the bottom panels of Fig.[3J). Each 
O V star was corrected for reddening according to the cor- 
responding cell E(B — V) and Ry values in Belikov et al. 
( 1999). From these stars we derive an absolute average dis- 
tance modulus of (m — M)o = 11.39 ± 0.46, correspond- 
ing to a distance from the Sun d Q = 1.92 ± 0.45 kpc. We 
use this value to set the 4 Myr Padova isochrone in the 
reddening-corrected, field-star decontaminated CMD of 
the central 6' of NGC 6611 in Fig.JHl Within uncertainties, 
the least- massive stars already in the MS of NGC 6611 
have mass in the range 3-5 M . 

Alternatively, we estimate the distance of NGC 6611 
by means of the MS entry point of the PMS stars. To do 
this we use the turn-on locus in the reddening-corrected, 
field-star decontaminated CMD of the central 6' in Fig. 03 
The PMS stars in Fig.OOturn out to be basically contained 
within the 200 Kyr and 4 Myr isochrones (Siess, Dufour, & 
Forestini 2000J), after applying an absolute distance mod- 
ulus of (m-M)o = 11.09 ± 0.30. This corresponds to 
a distance from the Sun of d© = 1.65 ± 0.27 kpc which, 
within the uncertainties, agrees with the previous value 
based on the OV stars. For completeness we also show 
in Fig. |5] evolutionary tracks up to the zero-age main se- 
quence (ZAMS) for PMS stars with masses 1, 2, 3 and 
5Mq. Similarly to the value based on the 4 Myr Padova 
isochrone, the low-mass MS of NGC 6611 implied by the 
PMS occurs in the range 3-5 M© . 

Based on the above estimates we adopt the average 
value d© = 1.8 ± 0.5 kpc as the distance of NGC 6611 
from the Sun. With this value the Galactocentric distance 
of NGC 6611 turns out to be dec = 6.3 ± 0.5 kpc, using 
8.0 kpc as the distance of the Sun to the center of the 
Galaxy (Reid I1993[) . The present value of the distance 
from the Sun agrees with that derived by Tadross et al. 
lj2UU2|l and is comparable to that of WEBDA (1.75 kpc), 
but falls ~ 0.5 kpc short with respect to the average value 
found in the literature. 

4. Cluster structure 

Structural parameters of NGC 6611 were derived by means 
of the radial density profile, defined as the projected num- 
ber of stars per area around the cluster center. In the case 
of this low-latitude cluster the contamination of the CMD 
by field stars must be taken into account in order to de- 
rive the intrinsic radial distribution of stars. We do this by 
considering the colour-magnitude filter shown in the inset 
of Fig. [5] which, within uncertainties, describes the cluster 
CMD morphology from the upper MS to the PMS stars. 
Before counting stars we applied the colour-magnitude 
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Field-star subtracted (0<R(')<6) 

Padova 4Myr: (m-M) =1 1 .39±0.45 

SiesslMyr: (m-M) =1 1.09+0.30 

Siess 4Myr: (m-M) =1 1 .09±0.30 
Siess 200Kyr: (m-M) =1 1 .0910.30 




100 - 



0.3 0.5 
(J-H) 

Fig. 5. The 4Myr Padova isochrone set to the reddening- 
corrected, field-star decontaminated CMD of the central 
6' of NGC6611 according to d Q = 1.92 ± 0.45 kpc. PMS 
model isochrones of Siess, Dufour, & Forestini 120001) with 
ages of 200 Kyr, 1 Myr and 4Myr are fitted to the PMS 
population resulting in d Q = 1.65±0.27kpc. Evolutionary 
tracks for PMS stars of masses 1, 2, 3and 5M Q are shown 
as dotted lines. 



filter to the differential-reddening corrected CMD of the 
cluster (including stars with distance to the center from 
0' to 23'), to discard stars with discrepant colours. The 
colour-magnitude filtering procedure has been previously 
applied in the analysis of the open clusters M 67 (Bonatto 
& BicaEOOS, N GC 188 (Bonatto, Bica & Santos Jr. lMHjl 
and NGC3680 (Bonatto, Bica & Pavani HUH^ . The ra- 
dial density profile was obtained by counting stars in- 
side concentric annuli with a step of 0.5' in radius. In 
the present case the background contribution level corre- 
sponds to the average number of stars included in the ring 
located at 17' < R < 22' (8.8 < R(pc) < 11.4), resulting 
ino-bg = 2.14±0.06 stars (')~ 2 = 7.8±0.2 stars pc" 2 . Note 
that the colour-magnitude filter encompasses both the MS 
and PMS stars. 

The resulting radial density profile of the MS + PMS 
stars is shown in Fig. El For absolute comparison between 
clusters we scale the radius in the abscissa in parsecs, and 




10.0 



Fig. 6. Radial density profile of the cluster NGC6611. 
Dashed line: 2-parameter King profile. Shaded region: stel- 
lar background level (average number of stars in the ring 
17' < R < 22'). The dimension of the core is indicated. 
The colour-magnitude filter used to discard stars with dis- 
crepant colours is shown in the inset, superimposed on the 
reddening-corrected and field-star decontaminated CMD 
of the central 6'. 



the number density of stars in the ordinate in stars pc~ 2 
using the distance derived in Sect. |31 The statistical sig- 
nificance of the profile is reflected in the ler Poisson error 
bars. 

Structural parameters of NGC6611 were derived by 
fitting the two-parameter King l|1966a|l surface density 
profile to the background-subtracted radial distribution 
of stars. The two-parameter King model essentially de- 
scribes the intermediate and central regions of normal 
clusters (King 1966b; Trager, King & Djorgovski 1995). 
The fit was performed using a nonlinear least-squares fit 
routine which uses the errors as weights. The best-fit so- 
lution is shown in Fig. superimposed on the observed 
radial density profile. Parameters derived are the King 
background-subtracted central density of stars ctok = 
35 ±6 stars (') -2 = 128 ±20 stars pc~ 2 , and the core radius 
Rcoro = 1-35 ± 0.15' = 0.7 ± 0.1 pc. Considering the radial 
density profile fluctuations with respect to the background 
level, we can define a limiting radius (Rii m ) for the cluster, 
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in the sense that for regions beyond Rn m , the null-contrast 
between cluster and background star density would pro- 
duce large Poisson errors and consequently, meaningless 
results. Thus, for practical purposes, the bulk of the clus- 
ter stars are contained within Rn m . For NGC6611 we es- 
timate a limiting radius Rii m = 12.5 ± 1.0' = 6.5 ± 0.5 pc. 

Despite the young age of NGC6611 the King profile 
provides a good fit of the stellar content in this cluster, 
within uncertainties. Since it follows from an isothermal 
(virialized) sphere, the close similarity of the radial dis- 
tribution of stars in NGC6611 with a King profile may 
suggest that the internal structure of this cluster (partic- 
ularly the core) has reached some level of energy equipar- 
tion. However, part of this effect appears to be related to 
formation processes as well (Sects. l5~2l and 

5. Luminosity and mass functions of the MS 

The rather populous nature of NGC 6611 provides an op- 
portunity to study the spatial distribution of LFs and MFs 
{(j)(m) = 4^) of the MS stars in such a young open clus- 
ter. 

Based on the King profile fit of the MS radial density 
profile of NGC 6611 (Fig.® we decided to study the MF 
in the following regions: (i) 0.0 < R(pc) < 0.7 (core), (ii) 
0.7 < R(pc) < 6.5 (halo), and (Hi) 0.0 < R(pc) < 6.5 
(overall). In order to maximize the significance of back- 
ground counts, we consider as offset field the outermost 
ring at 8.8 < R(pc) < 11.4, which lies ~ 2.3 pc beyond 
the limiting radius. 

Similarly to the structural analysis (Sect. the first 
step in the present analysis is to take into account the 
field-star contamination. In this case we apply the colour- 
magnitude filter (Fig. |5J) restricted to the MS stars (i.e. 
J < 11.4). The filtering process takes into account most 
of the field stars, leaving only a residual contamina- 
tion. We deal with this residual contamination statisti- 
cally by building LFs for each cluster region and offset 
field (17' < R < 22') separately. We consider the three 
2MASS bands independently when building the LFs. Note 
that the faint-magnitude limit of the MS is significantly 
brighter than that of the 99.9% Point Source Catalogue 
Completeness Limit (Sect.|3J|. We take J = 5 as the upper- 
MS limit (Fig.©. We build the LF of each 2MASS band 
by counting stars in magnitude bins from the respective 
faint magnitude limit to the upper MS, both for each 
cluster region and offset field. Considering that the solid 
angle of the offset field may be different from that of a 
given cluster region, we multiply the offset field LF by a 
numerical factor so that the solid angles match. The in- 
trinsic LF of each cluster region is obtained by subtract- 
ing the respective (i.e. solid angle-corrected) offset field 
LF from that of the cluster region. Finally, the intrinsic 
LFs are transformed into MFs using the mass-luminosity 
relation obtained from the 4Myr Padova isochrone (the 
youngest one available) and the absolute distance mod- 
ulus (m — M)o = 11.39 (Sect. OJ- These procedures are 
repeated independently for the three 2MASS bands. The 




4 5 6 7 8 9 10 15 20 30 



m(M D ) 

Fig. 7. Mass functions of the MS stars (empty circles) in 
different spatial regions of NGC 6611. Each panel contains 
MFs derived from the J, H and K s 2MASS photometry. 
MF fits {(j)(m) oc m ( 1+ x)) are shown as dashed lines, and 
the respective MF slopes are given. The insets show the 
MS + PMS field-star subtracted LFs in each spatial region. 



final MF of a given cluster region is produced by combin- 
ing the J, H and Ks MFs into a single MF. The resulting 
core, halo and overall MFs of NGC 6611, covering the mass 
range 5 < m(M ) < 25, are shown in Fig.[7| To these MFs 
we fit the function (p(m) oc m~( 1+x \ The resulting fits are 
shown in Fig. [7| and the MF slopes are given in col. 4 of 
Table [U 

We show in each panel of Fig. the corresponding 
field-star subtracted LFs, where the relative contribution 
of the PMS (Mj > -0.12) with respect to the MS can be 
evaluated. 

Both the overall and halo MFs are slightly steeper than 
a typical Salpeter 1|1955|> IMF (x — 1.35); however, the 
core MF turned out to be flat, with \ = 0.62 ± 0.16. The 
difference in MF slope between the core and halo (x — 
1.52 ± 0.13, panel (b) of Fig. |7|), may be a consequence of 
mass segregation in the core. We discuss this point further 
in Sect.Q 
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5.1. Total observed mass 

We provide in Table parameters derived from the MFs 
of each spatial region considered. The number of PMS 
stars in the core, halo and overall regions was estimated 
upon integration of the respective field star-subtracted 
LFs (Fig. Q) for Mj > -0.12. Considering the observed 
distribution of PMS stars and the respective evolutionary 
tracks (Fig. |5J) and the fact that MFs in general increase 
in number for the subsolar-mass range, for simplicity we 
assume an average PMS mass of 1 M Q . The resulting num- 
ber and mass of PMS stars are in cols. 2 and 3, respectively 
The observed number of MS stars and corresponding mass 
(cols. 5 and 6, respectively) are derived by integrating the 
MF through the mass range 5-25 M . We add to these 
the corresponding values of the number and mass of PMS 
stars to derive the total number of observed stars (col. 7), 
observed mass (col. 8), projected mass density (col. 9) and 
mass density (col. 10). 

To complete the census of the mass stored in observed 
stars in NGC 6611 we must include as well the stars more 
massive than 25 M , which although present in the pho- 
tometry, have not been counted in the mass estimate ow- 
ing to the adopted MF upper cutoff (25 M Q ). Hillenbrand 
et al. (1993) give the number of stars in selected mass 
bins for massive stars: 5 stars with mass in the range 
25 < m(M Q ) < 40, 4 with 40 < m(M ) < 60, and 1 with 
60 < m(M ) < 85. Accordingly, the mass stored in stars 
more massive than 25 M amounts to f» 435 M , about 
50% of the mass implied by the MF of the main sequence 
range considered above (col. 6 of Table . For simplic- 
ity we add their mass to the overall region. To this effect 
we include an additional entry in Tabled to account for 
the massive stars, resulting in iums+pms ~ 1630 M . This 
value is ~ 4 times larger than the estimate of Tadross et al. 

The total observed mass of NGC 6611 is a factor 
w 2 larger than the mass estimated for the Trapezium 
cluster (Lada & Lada I2003[) . It should be noted that 
the above mass and density estimates (cols. (8) - (10) in 
Table are in fact lower limits, since we are not taking 
into account stars less massive than 5M , a mass range 
for which we have no information on the behaviour of the 
MF slope. 

5.2. Dynamical state of NGC 6611 

The overall MF slope (\ — 1.45 ± 0.12) in the mass range 
5 < m(M ) < 25 is similar to that of a standard Salpeter 
(x — 1.35) IMF. Besides, the MF slope presents large 
variations in the inner regions, being flat (x = 0.62±0.16) 
in the core and rather steep (x = 1-52 ± 0.13) in the halo 
(Tableland Fig.EJ). In older clusters this fact reflects the 
effects of large-scale mass segregation, in the sense that 
low-mass stars originally in the core are transferred to the 
cluster's outskirts while massive stars accumulate in the 
core. This produces a flat MF in the core and a steep one 
in the halo (see, e.g. Bonatto & Bica 2005). 



Mass segregation in a star cluster scales with the 
relaxation time, defined as t ro i ax = gl ^ N t cross , where 
t C ross = R/cv is the crossing time, N is the (total) num- 
ber of stars and cr v is the velocity dispersion (Binney & 
Tremaine ll987fl . The characteristic time scale in which a 
cluster reaches some level of kinetic energy equipartition 
is described by t re i ax , when massive stars sink to the core 
and low-mass stars are transferred to the halo. Assuming 
a typical a v » 3kms _1 (Binney & Merrifield ll998|l we ob- 
tain for the overall cluster t re iax ~ 20 ± 5 Myr, and for the 
core t re ia X ~ 0.7 ± 0.2 Myr. Considering that NGC 6611 
is ~ 1.3 Myr old (Sect. [SJ, the age of this cluster is ~ 2 
times larger than t re i a x(core). Thus, the presence of mass 
segregation and consequently some degree of MF slope 
flattening can be expected in the core. However, the ratio 
cluster age to t re i ax drops to ~ 0.07 for the overall cluster, 
which is consistent with the typical Salpeter slope and the 
absence of important mass segregation. 

Although OB stars in NGC 6611 host more compan- 
ions (for binaries with mass ratios q > 0.1) than solar- type 
field stars (Duchene et al. l2001J) . the expected influence of 
binaries in the MF would be to steepen the slope of a 
standard Salpeter IMF to x ~ 1-7, and not to flatten it 
out (Sagar & Richtler 115511 Kroupa lMHjl . 

6. Comparison with clusters in different 
dynamical-evolution states 

In the previous sections we derived a series of parame- 
ters related to the structure and dynamical evolution of 
NGC 6611. At this point it may be useful to check how 
this young cluster fits in the context of open clusters with 
different ages and in more advanced dynamical states. 
Bonatto & Bica ( 2005 ) presented a systematic analysis of 
a set of open clusters with ages in the range 70- 7 000 Myr 
and masses from 400 - 5 300 M . The methodology used 
follows the same lines as that used the present paper. As 
a result, a set of uniform parameters related to the struc- 
ture (core and overall radii, mass and density), dynamical 
state (core and overall MF slopes, evolutionary parameter 
r = age/t re iax), as well as age and Galactocentric distance 
of open clusters was obtained. Some correlations among 
these parameters were verified, and a separation of mas- 
sive (m > 1000M Q ) and less-massive (m < 1000M o ) 
clusters was observed. The core and overall dynamical- 
evolutionary parameter of NGC 6611 is given in col. 11 of 
Tabled 

In Fig.|Slwe compare NGC 6611 with the set of nearby, 
well-studied older clusters (Bonatto & Bica 2005) in terms 
of structural parameters. Note that the mass and den- 
sity values of NGC 6611 are lower limits (Sect. I5.1|l . The 
limiting radius of NGC 6611 fits in the low- limit correla- 
tion of Rum with age (panel (d)). Although with a larger 
scatter, the same is observed in the correlation of R core 
with age (panel (g)). The core density of NGC 6611 fol- 
lows the trend presented by massive clusters for young 
ages (panel (a)). A similar trend is seen for the overall 
density (panel (h)). In Galactocentric distance the limit- 
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Table 1. Parameters derived from the MFs and PMS stars. 





PMS 


Observed MS 






Observed MS + PMS 




T 


Region 


N* 


m 


X5-25 


N* 


m bs 


N* 


niMS+PMS 




P 




(pc) 


(10 2 stars) 


(1O 2 M ) 




(10 2 stars) 


(1O 2 M ) 


(10 2 stars) 


(1O 2 M ) 


(M pc- 2 ) 


(M Q pc" 3 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Core 


1.0 ±0.3 


1.0 ±0.3 


0.62 ±0.16 


0.3 ±0.2 


2.3 ± 1.3 


1.2 ±0.4 


4.1 ± 1.6 


265 ± 107 


284 ± 115 


1.8 ±0.6 


Halo 


2.6 ±0.8 


2.6 ±0.8 


1.52 ± 0.13 


0.8 ±0.3 


6.9 ±2.9 


3.4 ±0.9 


9.6 ±3.1 


7.3 ±2.3 


0.83 ± 0.27 




Overall 


3.6 ± 1.2 


3.6 ± 1.2 


1.45 ±0.13 


0.9 ±0.3 


8.4 ±3.2 


4.5 ± 1.2 


12.0 ±3.4 


9.0 ±2.6 


1.0 ±0.3 


0.067 ±0.022 


Overall^ 


3.6 ± 1.2 


3.6 ± 1.2 


1.45 ±0.13 


1.0 ±0.3 


12.7 ±3.5 


4.6 ± 1.3 


16.3 ±3.5 


12.3 ± 2.6 


1.4 ±0.3 


0.066 ±0.022 



Table Notes. Observed mass range: 5 — 25 M . Col. 4 gives the MF slope derived for the MS stars. Col. 11: dynamical-evolution 
parameter r = age/t rc i ax . (f): includes stars more massive than 25 M . 
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Fig. 8. Relations involving structural parameters of 
open clusters. Filled circles: clusters more massive than 
1000M©. Open circles, m < 1000M Q . Filled triangle: 
NGC6611. Dashed lines: least-squares fits to the nearby 
clusters (see text). Shaded areas: la borders of the least- 
squares fits. Arrows indicate lower-limit estimates of mass 
and density for NGC6611. 



ing radius of NGC 6611 helps defining a correlation (panel 
(b)) in the sense that clusters at larger dec tend to be 



larger, which agrees with the results of Lynga (1982) and 
Tadross et al. l)2002[l . The large scatter in panel (e) pre- 
cludes any conclusion with respect to a dependence of the 
core radius with dec- NGC 6611 fits well in the tight cor- 
relations of core and overall mass (panel (i) ; least-squares 
fit: M coro = (14.71 ± 10.30) + (0.08 ± 0.01)M ovorall with a 
correlation coefficient CC — 0.92), and core and limiting 
radii (panel (f); R Um = (1.05 ± 0.45) ± (7.73 ± 0.66)R coro 
with CC = 0.95). Finally, in Bonatto, & Bica (2005} we 
observed that massive and less-massive clusters follow par- 
allel, different paths in the plane core radius x overall mass. 
NGC 6611 consistently fits in the the massive clusters path 
(panel (c)). 



NGC 6611 is compared to older clusters in terms 
of dynamical-evolution related parameters in Fig. [§J 
NGC 6611 fits in the tight correlations of core and overall 
MF slopes (panel (d); Xcorc = (-2.07 ± 0.46) + (1.74 ± 
0.35) x Xovcraii, with CC = 0.87), and core and over- 
all evolutionary parameters (panel (c); r coro = (63.3 ± 
9.6) x T OV erall) with CC = 0.91). The young age and small 
Galactocentric distance of NGC 6611 suggests a corre- 
lation of dec with age (panel (g)), which would agree 
with the result of Lynga l|1982fl . However, we note that 
this trend is mostly based on a single point (NGC 6611) 
and thus, a larger cluster sample within the central sec- 
tor is necessary to establish the statistical signficance 
of the correlation. The Salpeter-like overall MF slope of 
NC6611 is consistent with the overall mass and age re- 
lation when compared to other massive clusters (panels 
(a) and (i), respectively). The scatter in panel (h) does 
not allow inferences on the relation of core MF slope 
with mass. Finally, with respect to the evolutionary pa- 
rameters of NGC 6611, the overall MF slope (panel (b); 

Xovcraii = (1-5 ± 0.1) - (1.9 ± 0.6) exp - ( 16 T ^ 3 1 - 1 ) , with 
CC = 0.89) follows the relations suggested by massive 
clusters with small values of r. However, the core MF slope 
is flatter than that expected for its value of r (panel (e); 
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Fig. 9. Relations involving evolutionary parameters of 
open clusters. Symbols as in Fig. |S] 

Xcorc = (1.17 ± 0.23) - (3.0 ± 0.7) exp - , with 

CC = 0.82) and age (panel (f)). 

We conclude that both structurally and dynamically, 
the overall parameters of NGC6611 fit in the relations 
defined by older, more-massive open clusters, in the loci 
expected of a very young cluster. 

The core of NGC 6611 presents a MF significantly flat- 
ter than that in the halo, a situation typical of older clus- 
ters with important mass segregation. Part of this effect 
can be accounted for by the fact that the core relaxation 
time is about half the cluster age. 

7. Concluding remarks 

In this paper we analyzed the structure and mass distri- 
bution of MS and PMS stars in the ionizing open clus- 
ter NGC 6611. The analysis was based mostly on J, H 
and Ks 2MASS photometry corrected for the variable 
reddening across the region, from which we built radial 
density profiles and mass functions. The average colour 
excess in a region reaching R « 23' from the center of 
NGC 6611 resulted in E(J - H) = 0.25, corresponding to 
E(B — V) = 0.8. To uncover the intrinsic CMD morphol- 
ogy we applied a field-star decontamination procedure to 



the reddening-corrected data. The resulting CMDs present 
a MS restricted to stars more massive than « 5 Mq and a 
collection of PMS stars. 

Based on the fraction of stars with Kg excess we es- 
timated an age of - 1.3 ± 0.3 Myr for NGC 6611. We de- 
rived a distance from the Sun of d© = 1.8 ± 0.5 kpc, con- 
sidering the values calculated from OV stars with previ- 
ous spectroscopic data and PMS tracks. Our value of d Q 
agrees with those given by Tadross 1)2002(1 and WEBDA. 
However, it locates NGC 6611 ~ 0.5 kpc closer to the Sun 
than the average value given in previous works. The di- 
rect method of distance determination using stars is fun- 
damental for the definition of the Galactic structure. Note 
that available kinematic distances of the related Hn re- 
gion Sh2-49 are d Q = 2.95 kpc (Georgelin, & Georgelin 
IHJ70}) and d Q = 2.2 kpc (Blitz, Fich, & Stark Egg) , thus 
larger than the present value. 

King model fits to the radial density profiles including 
MS and PMS stars produced a core radius R CO rc = 0.70 ± 
0.08 pc with a limiting radius Rn m » 6.5 ± 0.5 pc. The 
projected density of stars in the cluster center is ~ 128 ± 
20 stars pc -2 . 

Considering the field star-subtracted LFs, we esti- 
mated that the number of PMS stars in NGC 6611 is 
~ 362. The total observed mass locked up in PMS and 
MS stars amounts to hims+pms = (1-6 ± 0.3) x 10 3 M Q . 
This is a lower-limit estimate since we are not taking into 
account stars less massive than 5 M Q . 

The core MF of the MS stars in NGC 6611 is flat with 
a slope x — 0-62 ± 0.16, while in the halo it steepens to 
X = 1.52 ± 0.13. The overall MF of NGC 6611 is similar 
to a Salpeter IMF with x = 1-45 ± 0.12. 

Compared to older open clusters in terms of structural 
and dynamical-cvolution-related parameters, the overall 
cluster behaves as a massive open cluster at young ages. 

The spatial variation of MF slopes, being flat in the 
core and steeper in the halo, implies that mass segrega- 
tion has already affected the mass distribution in the core, 
since the relaxation time t re i a x(core) ~ 0.73 Myr corre- 
sponds to ~ 50% of the cluster age. For the whole cluster 
treiax (overall) ~ 20 Myr ~ 15 x cluster age, which means 
that mass segregation did not have time to redistribute 
stars in large scale throughout the body of such a young 
cluster. Accordingly, the flattening degree in the core MF 
of NGC 6611 seems to be related as well to initial condi- 
tions, probably associated to the fragmentation of the par- 
ent molecular cloud itself, with more massive proto-stars 
preferentially located in the central parts of the cloud. 
This agrees with Kroupa (2002J with respect to the ob- 
served degree of mass segregation in clusters younger than 
a few Myr. 
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